Prognostic impact of a lymphocyte activation-associated gene signature in GBM based on transcriptome analysis.

BACKGROUND: Glioblastoma multiforme (GBM) is a highly, malignant tumor of the primary central nervous system. Patients diagnosed with this type of tumor have a poor prognosis. Lymphocyte activation plays important roles in the development of cancers and its therapeutic treatments. OBJECTIVE: We sought to identify an efficient lymphocyte activation-associated gene signature that could predict the progression and prognosis of GBM. METHODS: We used univariate Cox proportional hazards regression and stepwise regression algorithm to develop a lymphocyte activation-associated gene signature in the training dataset (TCGA, n = 525). Then, the signature was validated in two datasets, including GSE16011 (n = 150) and GSE13041 (n = 191) using the Kaplan Meier method. Univariate and multivariate Cox proportional hazards regression models were used to adjust for clinicopathological factors. RESULTS: We identified a lymphocyte activation-associated gene signature (TCF3, IGFBP2, TYRO3 and NOD2) in the training dataset and classified the patients into high-risk and low-risk groups with significant differences in overall survival (median survival 15.33 months vs 12.57 months, HR = 1.55, 95% CI [1.28-1.87], log-rank test P < 0.001). This signature showed similar prognostic values in the other two datasets. Further, univariate and multivariate Cox proportional hazards regression models analysis indicated that the signature was an independent prognostic factor for GBM patients. Moreover, we determined that there were differences in lymphocyte activity between the high- and low-risk groups of GBM patients among all datasets. Furthermore, the lymphocyte activation-associated gene signature could significantly predict the survival of patients with certain features, including IDH-wildtype patients and patients undergoing radiotherapy. In addition, the signature may also improve the prognostic power of age. CONCLUSIONS: In summary, our results suggested that the lymphocyte activation-associated gene signature is a promising factor for the survival of patients, which is helpful for the prognosis of GBM patients.

Transcriptomic heterogeneity of driver gene mutations reveals novel mutual exclusivity and improves exploration of functional associations.

BACKGROUND: Lung adenocarcinoma (LUAD), as the most common subtype of lung cancer, is the leading cause of cancer deaths in the world. The accumulation of driver gene mutations enables cancer cells to gradually acquire growth advantage. Therefore, it is important to understand the functions and interactions of driver gene mutations in cancer progression. METHODS: We obtained gene mutation data and gene expression profile of 506 LUAD tumors from The Cancer Genome Atlas (TCGA). The subtypes of tumors with driver gene mutations were identified by consensus cluster analysis. RESULTS: We found 21 significantly mutually exclusive pairs consisting of 20 genes among 506 LUAD patients. Because of the increased transcriptomic heterogeneity of mutations, we identified subtypes among tumors with non-silent mutations in driver genes. There were 494 mutually exclusive pairs found among driver gene mutations within different subtypes. Furthermore, we identified functions of mutually exclusive pairs based on the hypothesis of functional redundancy of mutual exclusivity. These mutually exclusive pairs were significantly enriched in nuclear division and humoral immune response, which played crucial roles in cancer initiation and progression. We also found 79 mutually exclusive triples among subtypes of tumors with driver gene mutations, which were key roles in cell motility and cellular chemical homeostasis. In addition, two mutually exclusive triples and one mutually exclusive triple were associated with the overall survival and disease-specific survival of LUAD patients, respectively. CONCLUSIONS: We revealed novel mutual exclusivity and generated a comprehensive functional landscape of driver gene mutations, which could offer a new perspective to understand the mechanisms of cancer development and identify potential biomarkers for LUAD therapy.

Systematic Assessment of Transcriptomic Biomarkers for Immune Checkpoint Blockade Response in Cancer Immunotherapy.

BACKGROUND: Immune checkpoint blockade (ICB) therapy has yielded successful clinical responses in treatment of a minority of patients in certain cancer types. Substantial efforts were made to establish biomarkers for predicting responsiveness to ICB. However, the systematic assessment of these ICB response biomarkers remains insufficient. METHODS: We collected 22 transcriptome-based biomarkers for ICB response and constructed multiple benchmark datasets to evaluate the associations with clinical response, predictive performance, and clinical efficacy of them in pre-treatment patients with distinct ICB agents in diverse cancers. RESULTS: Overall, "Immune-checkpoint molecule" biomarkers PD-L1, PD-L2, CTLA-4 and IMPRES and the "Effector molecule" biomarker CYT showed significant associations with ICB response and clinical outcomes. These immune-checkpoint biomarkers and another immune effector IFN-gamma presented predictive ability in melanoma, urothelial cancer (UC) and clear cell renal-cell cancer (ccRCC). In non-small cell lung cancer (NSCLC), only PD-L2 and CTLA-4 showed preferable correlation with clinical response. Under different ICB therapies, the top-performing biomarkers were usually mutually exclusive in patients with anti-PD-1 and anti-CTLA-4 therapy, and most of biomarkers presented outstanding predictive power in patients with combined anti-PD-1 and anti-CTLA-4 therapy. CONCLUSIONS: Our results show these biomarkers had different performance in predicting ICB response across distinct ICB agents in diverse cancers.

Identification of a Six-lncRNA Signature With Prognostic Value for Breast Cancer Patients.

Breast cancer (BRCA) is the most common cancer and a major cause of death in women. Long non-coding RNAs (lncRNAs) are emerging as key regulators and have been implicated in carcinogenesis and prognosis. In this study, we aimed to develop a lncRNA signature of BRCA patients to improve risk stratification. In the training cohort (GSE21653, n = 232), 17 lncRNAs were identified by univariate Cox proportional hazards regression, which were significantly associated with patients' survival. The least absolute shrinkage and selection operator-penalized Cox proportional hazards regression analysis was used to identify a six-lncRNA signature. According to the median of the signature risk score, patients were divided into a high-risk group and a low-risk group with significant disease-free survival differences in the training cohort. A similar phenomenon was observed in validation cohorts (GSE42568, n = 101; GSE20711, n = 87). The six-lncRNA signature remained as independent prognostic factors after adjusting for clinical factors in these two cohorts. Furthermore, this signature significantly predicted the survival of grade III patients and estrogen receptor-positive patients. Furthermore, in another cohort (GSE19615, n = 115), the low-risk patients that were treated with tamoxifen therapy had longer disease-free survival than those who underwent no therapy. Overall, the six-lncRNA signature can be a potential prognostic tool used to predict disease-free survival of patients and to predict the benefits of tamoxifen treatment in BRCA, which will be helpful in guiding individualized treatments for BRCA patients.

Single-cell RNA-seq reveals the invasive trajectory and molecular cascades underlying glioblastoma progression.

Glioblastoma (GBM) is the most common and aggressive primary brain tumor, in which GBM stem cells (GSCs) were identified to contribute to aggressive phenotypes and poor prognosis. Yet, how GSCs progress to invasive cells remains largely unexplored. Here, we revealed the cell subpopulations with distinct functional status and the existence of cells with high invasive potential within heterogeneous primary GBM tumors. We reconstructed a branched trajectory by pseudotemporal ordering of single tumor cells, in which the root showed GSC-like phenotype while the end displayed high invasive activity. Thus, we further determined a path along which GSCs gradually transformed to invasive cells, called the 'stem-to-invasion path'. Along this path, cells showed incremental expression of GBM invasion-associated signatures and diminishing expression of GBM stem cell markers. These findings were validated in an independent single-cell data set of GBM. Through analyzing the molecular cascades underlying the path, we identify crucial factors controlling the attainment of invasive potential of tumor cells, including transcription factors and long noncoding RNAs. Our work provides novel insights into GBM progression, especially the attainment of invasive potential in primary tumor cells, and supports the cancer stem cell model, with valuable implications for GBM therapy.

Identification of a five-lncRNA signature for predicting the risk of tumor recurrence in patients with breast cancer.

Long non-coding RNAs (lncRNAs) are a major class of non-coding RNAs, and the functional deregulations of lncRNAs have been shown to be associated with the development and progression of BC. In this work, we conduct an integrative analysis on five re-annotated lncRNA expression datasets from the Gene Expression Omnibus (GEO) which included a total of 891 BC samples. We identified a five-lncRNA signature that was significantly associated with DFS in the training cohort of 327 patients. We found the five-lncRNA signature could effectively stratify patients in the training dataset into high- and low-risk groups with significantly different DFS (p = 3.29 × 10-5 , log-rank test). The five-lncRNA signature was effectively validated in four independent cohorts, and prognostic analysis results showed that the five-lncRNA signature was independent of clinical prognostic factors, such as BC subtypes and adjuvant treatments. Furthermore, GSEA suggested that the five-lncRNA signature was involved in BC metastasis-related pathways. Our findings indicate that these five lncRNAs may be implicated in BC pathogenesis, and further, these lncRNAs may potentially serve as novel candidate biomarkers for the identification of BC patients at high risk for tumor recurrence.

LncNetP, a systematical lncRNA prioritization approach based on ceRNA and disease phenotype association assumptions.

Our knowledge of lncRNA is very limited and discovering novel disease-related long non-coding RNA (lncRNA) has been a major research challenge in cancer studies. In this work, we developed an LncRNA Network-based Prioritization approach, named "LncNetP" based on the competing endogenous RNA (ceRNA) and disease phenotype association assumptions. Through application to 11 cancer types with 3089 common lncRNA and miRNA samples from the Cancer Genome Atlas (TCGA), our approach yielded an average area under the ROC curve (AUC) of 83.87%, with the highest AUC (95.22%) for renal cell carcinoma, by the leave-one-out cross validation strategy. Moreover, we demonstrated the excellent performance of our approach by evaluating the influencing factors including disease phenotype associations, known disease lncRNAs and the numbers of cancer types. Comparisons with previous methods further suggested the integrative importance of our approach. Taking hepatocellular carcinoma (LIHC) as a case study, we predicted four candidate lncRNA genes, RHPN1-AS1, AC007389.1, LINC01116 and BMS1P20 that may serve as novel disease risk factors for disease diagnosis and prognosis. In summary, our lncRNA prioritization strategy can efficiently identify disease-related lncRNAs and help researchers better understand the important roles of lncRNAs in human cancers.